Photon generation through decoherence in cavity QED: analytical analysis 
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The Rabi Hamiltonian describes the interaction between a two-level atom and a single mode of the 
quantized electromagnetic field. When the system is subject to the Markovian atomic dephasing 
reservoir, the anti-rotating term leads to the photon generation from vacuum. In the presence 
of Markovian damping reservoirs, the asymptotic mean photon number is higher than the thermal 
photon number expected in the absence of the anti-rotating term. We obtain approximate analytical 
expressions in the asymptotic regime for the photon creation rate in the pure dephasing case and 
the mean photon number and the atomic population inversion in the general case. Our analytical 
results are valid in the small mean photon number limit and they were tested by numerical analyses. 
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The Rabi Hamiltonian (RH) [l[ is the simplest and the 
most used Hamiltonian deduced from first principles to 
describe the interaction between a two-level atom and a 
single mode of the quantized electromagnetic (EM) field 
'2). It reads (h = 1) 

UJ Q 



H = con 



V2gxa x , 



(1) 



where u> and coq are the cavity and atomic transition fre- 
quencies, respectively, and g is the coupling constant. 
The atomic operators are 



where a a 



(e\ - Is) (g\ 



(g\ and er_ = a J + , with \g) and |e) de- 

The 



noting the ground and excited states, respectively, 
cavity field quadrature operators are 



V2 



P = 



V2i 



where a (od) is the annihilation (creation) operator of the 
EM field and n — a^a is the photon number operator. 

The exact solution of the RH is lacking to the present 
day, so usually one performs the RWA by which the 

so called anti-rotating term g (ader_|_ + aer_) is neglected. 
The resulting approximate Hamiltonian is known as the 
Jaynes-Cummings Hamiltonian (JCH) 0] and has an el- 
egant exact solution 0, H, [H, @ , which led to a predic- 
tion of a variety of purely quantum phenomena, such 
as collapse and revival of the atomic inversion Q, Rabi 
oscillations Q, squeezing @, non-classical states, such 
as the Schrodinger cat-like state [l(| and Fock states 
, and the atom-atom or atom-field entanglement [12] • 
Many of these phenomena were observed experimentally 
in the last decades, thereby giving an experimental val- 
idation test of the RWA. Nowadays, JCH is the main 
analytic tool to analyse cavity quantum electrodynamics 
(QED) 0, [l3j], including protocols for quantum informa- 
tion processing in systems composed by an effective two- 
level atom (qubit) and photons or phonons. As exam- 
ples, one can cite the already implemented architectures, 



such as circuit qedEI EE El, E3 , cavit y Q ED 03 El 
and trapped ions [19(, or novel proposals, such as polar 
molecules coupled to stripline resonators [20( , mechanical 
resonators coupled to an electronic spin qubit [2l| . etc. 

The range of validity of RWA has been studied for a 
long time and the common sense @, H, @ is that it is 
valid for a weak field amplitude, small g and small de- 
tuning |A| <C 1, where A = luq — uj. Several numeri- 
cal studies exemplified the deviation of the dynamics of 
the RH compared to the one expected from the JCH 
[S HI HI HlHl, thereby demonstrating the breakdown 
of RWA in specific regimes. Alternative approximations 
have been also su gge sted in order to increase the validity 
of the RWA [H, H3|. Moreover, it was shown recently 
that the antirotating term is responsible for the photon 
generation from vacuum due to an analog of the Dynami- 
cal Casimir Effect [28| in non-stationary cavity QED sys- 
tems m, us m ■ 

The majority of the previous studies on the role of the 
anti-rotating term was performed in the idealized closed 
system approach, when the atom-cavity system is iso- 
lated from the environment. Although the dynamics of 
the JCH has been also extensively studied in the pres- 
ence of dissipative environments [15|, l32|, |33|, [3J| , the RH 
has not received much attention in open system dynam- 
ics until recently [Uj]. However, it seems that novel and 
unexpected phenomena appear when one combines the 
anti-rotating term with the dissipation induced by sev- 
eral kinds of environments [35( . 

Recently, we studied numerically the dynamics of the 
cavity QED system (35|, described by the RH, subjected 
to the action of damping and dephasing reservoirs. Un- 
der the standard Born-Markovian and the weak system- 
reservoir coupling approximations [32l l3ol \3l\ , the system 
dynamics is governed by the master equation 



dp 
di 



= -i [H,p]+C(p) 



(2) 



where H is the RH ([T]) and the dissipation superoperator 



2 



C (p) is given by 

C (p) = C a (p) + C f (p) + C da (p) + C df (p) , (3) 
with the standard definitions [36| 

L a Go) = 7 (n t + 1)V[j-] + jn t V[a+] 
Cf(p) = n(n t + l)V[a] + Kn t V[a ] ] 



where 



£da (p) 



l P h 



2 1 1 



Above we used the short notation [15| for the Lindblad 
superoperator V[L]p = (2Lptf - L^Lp - pL^L) /2. The 
superoperators C a (p) and C / (p) describe the effects of 
the thermal reservoirs (with mean photon number nt) on 
the atom and the field, respectively, where 7 (k) is the 
atom (cavity) relaxation rate. Another source of deco- 
herence are the phase damping reservoirs acting on the 
atom (field), represented by Cda (p) (Cdf (p)), where j ph 
(T p h) is the atomic (cavity) pure dephasing rate. Usu- 
ally, the cavity dephasing is small compared to other dis- 
sipative channels in circuit QED, so it is neglected [38| . 
However, due to the measurement back-action (39|, the 
effective cavity dephasing rate can become large in mi- 
crowave cavity QED, in which the field is continuously 
measured via quantum non-demolition photon counting 
using non-resonant Rydberg atoms [ipj . 

In (35j we showed numerically that in the presence of 
only atomic dephasing reservoir (T p h = k = 7 = 0), 
there is asymptotically a linear photon growth as func- 
tion of time for any initial state |V , o)i even for atom 
and field being initially in their respective ground states, 
lV>o) = \g, 0). This occurs due to the combination of the 
atomic dephasing and the anti-rotating term. Moreover, 
for non-zero k and 7 the mean photon number attains a 
stationary value greater than the thermal photon number 
n t in the reservoir. The atomic population inversion also 
achieves a stationary value above the one expected from 
the JCH. Here we evaluate analytically the photon gen- 
eration rate for the pure dephasing case and estimate the 
asymptotic mean photon number and the atomic popu- 
lation inversion when all the sources of loss are present 
in the weak coupling regime and small photon number 
limit. 

From © the Heisenberg equations of motion for the 
mean photon number (n) and the atomic population in- 
version (a z ) read 



(n) = -V2g (pa x ) - K (n) + nn t (4) 
(a z ) = 2V2. 9 (xa v ) - 1 [1 + {2n t + 1) (a z )] . (5) 
The equations for the higher order dynamic variables are 



(pa^ = - (xa x ) - uj (pa y ) - x (pc*) - V2g 



xa x j = (pa x ) - uj a (xa y ) - x (xa x ) 
(p a y) = w o (p<y x ) - (xo y ) - x (p<7 v ) - V2ga (6) 
(xa y ^ = (pa y ) + uj (xa x ) - x (x<r y ) - V2g(, 



X = 7 P h + T ph + 77 + 7 [nt + - 



(7) 



The system of equations ([6]) is not closed because of the 
dynamic variables ( = (2x 2 a z ) and a = ((xp + px)a z ), 
which obey the corresponding differential equations. 
Therefore, this system of equations cannot be integrated 
exactly, although some numerical methods for its solution 
based on semi-Lie algebra have been proposed [24J. For 
simplicity, from now on we shall neglect the pure cavity 
dephasing, setting T p h = 0. 

First, we estimate the photon generation rate for the 
initial state \g, 0) in the absence of damping (k = 7 = 0), 
under the experimentally realistic weak coupling (g <C 1) 
and low temperature (n t <C 1) regimes. While the 
mean number of photons generated through decoherence 
is small, (n) <C 1, the probability of the state \g, 0) 
is high, so we assume C — (9, 0\^ x2(J z\g, 0) = —1 and 
a ~ (g, 0| (xp + px) a z \g, 0) = in order to make the sys- 
tem of equations © solvable. To test the validity of this 
assumption, we studied numerically the asymptotic val- 
ues C a an d OL a as function of A, g and j ph in the range 
of parameters we are interested in. From now on we set 
bj = 1. As shown in Fig. 1, £ a and a a are always close 
to —1 and 0, respectively, as expected. 

Thus, for (n) <C 1 the asymptotic values of the quan- 
tities appearing on the RHS of Eqs. ((4]) and ([5]) are 



(P Cr ^)a 



-y/2- 



91 ph 



A 



9 ' 

Tph 



where we defined A + = uj + luq. Substituting (pa x ) a into 
Eq. ((4]) we find that asymptotically the photon creation 
rate (h)„ attains a constant value 



27; 



ph 



(8) 



For 7 p/l <C A + , (n) a is proportional to j ph and g 2 , and 
inversely proportional to A^_ , as observed numerically in 
[35l ]. In Fig. 2 we compare the numerical values of (n) a 
to the formula {8} as function of j ph , g and A + . We 
see that the approximate formula scales correctly with 
the system parameters, but differs slightly from the nu- 
merical values due to the deviations of £ and a from — 1 
and 0, respectively, which were neglected in our approx- 
imation. Nevertheless, the simple formula © gives the 
correct order of magnitude of the photon generation rate 
through decoherence. 

Applying the same procedure to the general case with 
damping, in the limit (n) < 1 we obtain the following 
stationary values for the mean photon number (n) and 
the atomic population inversion (a z ) 
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TABLE I: Numerical values N and S and upper bounds 
(N > ,S > ) given by inequalities (|ll|l - (|12|l for different values 
of y ph , 7 and k for T ph = A = 0, g = 2 • 10~ 2 , n t = 0. The 
lower theoretical bounds are iV< = 10 -4 and S< = 4 • 10" 4 . 



where 



e = 



A 



We point out that in the special case uiq = one can 
easily get the exact expressions for (h) a and ("■) oc , which 
turn out to be precisely Eqs. J5]) and ([9]) with the equality 
sign, independently of the mean photon number. This is 
a consistency check of our treatment. 

We can also obtain the lower bounds for (n) and 
(<t z } 00 as follows. The minimum values of (n)^ and 
(cr z ) occur in the 'worst' scenario, when the density 
matrix is as close to \g, 0}(g, 0| as possible, so ( = —1 and 
a = hold almost exactly. Therefore, from ([9]) and (fl"0|) 
we obtain (n)^ > n*+6 and {<J z ) oa > — (2n t + 1) - +46 
and we get a simple inequality for the stationary values 



6 < (n)^ - th < 26 



46 < (a 2 



1 



2n t + 1 



< 46 



TK + 1/2)' 



(11) 
(12) 



From Eqs. (fTTj) and (fl"2")) we see that asymptoti- 
cally the anti-rotating term combined with dissipative 
losses creates at least 6 ~ (g/A + ) 2 photons above the 
thermal photon number n t . In the tables I and II we 
compare the numerical values of N = (n) — n t and 

S = ((Tz)^ + {2th + IT 1 to the lower (JV<, 5<) and up- 
per (JV>, 5>) bounds given by Eqs. (jlip and (p~2|) . show- 
ing that the lower bound is always satisfied. The upper 
bound holds in the majority of cases and is satisfied in the 
order of magnitude in all the simulations we performed. 



The expression (ITTT) could partially explain why the mean 
photon number observed in [40] is slightly higher than 
the expected thermal photon number - in that case, the 
cavity dephasing rate T p h induced by the measurement 
back-action can be significant, so one would expect the 
generation of photons from vacuum through decoherence. 

In summary, we studied analytically the phenomenon 
of photon generation through decoherence in the cav- 
ity QED system, described by the Rabi Hamiltonian, 

A+ 10 4 7V 10 4 7V<, 10 4 iV> 10* S 10 i S < , 10 i S > 
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TABLE II: The numerical values N and S, the lower (7V< , 5<) 
and the upper (iV> , 5> ) theoretical bounds given by inequal- 
ities (|ll|l - (|12p for different values of g and A+ for nt = and 
decay rates (F ph , -y ph , 7, «) = (0,0.1, 1, 1) • 10" 2 . 



coupled to Markovian dephasing and dissipative reser- 
voirs. We obtained a simple expression characterizing 
the asymptotic photon generation rate in the pure de- 
phasing case, which agrees with the numerical results 
previously obtained in Ref j35| . Moreover, we deduced 
approximate inequalities giving the lower and the up- 
per bounds for the asymptotic mean photon number and 
the atomic population inversion when the damping is 
present. These expressions were confirmed by numerical 
simulations, showing that the lower bound always holds 
and the upper bound is satisfied in the order of magni- 
tude. 6ur expressions also agree with the exact formula 
in the special case of null atomic frequency, providing one 
more consistency check of our treatment. We emphasize 
that our results are valid in the limit of small mean pho- 
ton number, which is precisely the situation one expects 
in realistic cavity QED systems. Therefore, our study 
demonstrates the importance of the anti-rotating term 
in the open system dynamics and gives an estimative of 
its influence on the experimentally observable quantities. 
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Figure Captions 

Fig. 1: Asymptotic values ( a and a a for T p h — 7 = 
k = as function of a) g for -f ph = 0.05 and A + = 2; b) 
-f ph for A + = 2 and g = 0.02; c) A + for g = 0.02 and 
7 pfc = 0.05. 

Fig. 2: Comparison between the numerical values of 
(n) a (dots) and Eq. ((5]) for r p /j = 7 = n = as function 
of a) g for 7ph = 0.05 and A+ = 2; b) lph for A+ = 2 
and g = 0.02; c) A+ for g = 0.02 and j h = 0.05. 
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